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Abstract 

We reconsider the results cocerning the extreme-quantum S = 1/2 square- 
lattice Heisenberg antiferromagnet with frustrating diagonal couplings (Ji — J2 
model) drawn from a comparison with exact-diagonalization data. A combined 
approach using also some intrinsic features of the self-consistent spin-wave theory 
leads to the conclusion that the theory strongly overestimates the stabilizing role 
of quantum flutcuations in respect to the Neel phase in the extreme-quantum case 
S = 1/2. On the other hand, the analysis implies that the Neel phase remains 
stable at least up to the limit J2/J1 = 0.49 which is pretty larger than some 
previous estimates. In addition, it is argued that the spin-wave ansatz predicts 
the existence of a finite range ( J2/J1 < 0.323 in the linear spin-wave theory ) 
where the Marshall-Peierls sigh rule survives the frustrations. 



1 Introduction 

The square-lattice Heisenberg antiferromagnet with antiferromagnetic next-nearest- 
neighbor couphngs (Ji — J2 model) produces a simple and, at the same time, an 
important example of a frustrated quantum spin system. The model is defined by 
the Hamiltonian 

i7 = JiJ^S.S. + JaES.S, , Ji,J2>0, (1) 

where the symbols < i,j > and [i,j] mean that the summations run over the 
nearest-neighbor and next-nearest-neighbor (diagonal) bonds, respectively. In 
what follows we put Ji = 1, a = J2/J1. 

Presently little is known about the ground-state properties of this model. In 
the classical limit S = 00 the Ji — J2 model has two phases: if a < 1/2 ,the ground 
state is a two-sublattice Neel state, whereas, if a > 1/2, the four-sublattice an- 
tiferromagnetic state is stable. At the classical transition point a = 1/2 the 
model is characterized by a great degree of classical degeneracy: all states with 
zero elementary-plaquette spins are energetically preferable. The quantum fluc- 
tuations, however, can drastically change this picture. In general, they are deter- 
mined by the microscopic structure of the model, and are expected to increase as 
S approaches the extreme-quantum limit S = 1/2, and/or the frustration becomes 
stronger. Already a simple linear spin- wave analysis reveals such a tendency [1]. 
In addition, the latter theory predicts the existence of a finite range around the 
classical phase boundary a = 1/2 where the classical long-range magnetic order 
is completely destroyed (for arbitrary S). However, the next-order terms in the 
large-S* expansion show logarithmic divergencies [2] connected to an additional 
softening of the spectrum at a = 1/2, thus making the first-order predictions, at 
least, questionable. This situation is characteristic for most of the studied frus- 
trated models. An open question is how to reconstruct the standard spin-wave 
expansion in order to avoid the mentioned difficulties. The Hartree-Fock type 
theories [2-4], which could in principle serve as a starting point for a systematic 
expansion, predict a first-order phase transition between the magnetically ordered 



phases without any intermediate phase. This picture is connected with the pre- 
dicted stabihzing role of quantum fluctuations in respect to the two-sublattice Neel 
order. Presently, however, it is not clear if these conclusions are characteristic, at 
least qualitatively, for the extreme-quantum system S* = 1/2 as well. 

Concerning the S = 1/2 case, at least two important issues, related to the 
ground-state phase diagram, remain unsettled: (i) the nature of the magnetically 
disordered phase, if any, in the strongly frustrated region; (ii) the location of the 
phase-transition boundary. The magnetically disordered spin-Peierls dimer state 
is preferable in a number of studies: 1) series expansions around dimer states 
[5], 2) 1/A^-expansion technique [6], 3) bond-operator techniques [7], 4) effective- 
action approaches leading to quantum nonlinear a-models [8], 5) numerical exact- 
diagonalization data [9,10]. However, each of the mentioned methods has its own 
defects, so that some other states (e.g., the chiral states [10,11,12]) seem to be 
possible candidates, as well. 

With regard to the location of phase boundary, here the estimates run in the 
large interval from a^ ~ 0.15 to ac ~ 0.6. The bound ac ~ 0.15 was obtained 
[13] by use of o"- model considerations combined with Schwinger-boson mean-field 
results for S = 1/2. On the other hand, the largest estimate ac ~ 0.6 is charac- 
teristic for the self-consistent theories [2-4]. Series of studies give values which are 
near the point ac = 0.4 [1,9,10]. 

The outlined ambiguity signals of a lack of reliable descriptions even in the 
weakly frustrated region where the two-sublattice Neel phase is expected to be sta- 
ble. Concerning the spin-wave theories, a way to test their quality gives the com- 
parison with numerical exact-diagonalization data. For the Ji — J2 model the first 
steps in this direction were made by Hirsch and Tang [14] based on Takahashi's 
idea [15] for a constrained spin- wave theory in low dimensions. These authors 
indicated that their theory systematically overestimates the effect of frustration 
in destroying the antiferromagnetic correlations (A^ = 10,16,26). Recently, Cec- 
catto, Gazza, and Trumper [16] have continued this line by use of Takahashi's 
self-consistent approach [17] adapted to the frustrated model [4]. A remarkable 



agreement with the exact results for a number of lattices (A^ = 10, 20, 26) was 
indicated, excluding, however, the most symmetrical lattice 4x4. 

In this paper we study the extreme-quantum system S = 1/2 and show that the 
existing exact-diagonalization results in combination with some intrinsic proper- 
ties of the self-consistent theory lead to the following conclusions: (i) The classical 
Neel state is stable at least up to the limit a* = 0.49. Notice that the estimate 
is quite larger than the result a ~ 0.4 mentioned above; (ii) The self-consistent 
spin-wave theory overestimates the role of quantum fluctuations in stabilizing the 
Neel state. This last conclusion also differ from previous considerations relying 
on a comparison with exact-diagonalization results for less symmetrical lattices 
A^ = 10, 20, 26 when the theory indeed gives exellent results. 

2 Comparison of the theory with exact-diagonahzation 
data 



2.1 Fitting to the exact-diagonalization results: A^ = 16 
lattice 

In reconsidering the previous exact-diagonalization data, it is easy to see a well- 
pronounced tendency, i.e., the self-consistent spin- wave theory gives good correla- 
tors for a number of less-symmetrical lattices (A^ = 10,20,26), whereas the most 
symmetrical 4x4 lattice is aside from this tendency. On the other hand, Hirsch- 
Tang's theory overestimates the effect of frustration for all lattices ( including 
the 4x4 lattice). In order to further check this observations, we present here new 
exact-diagonalization results for A^ = 18 and A^ = 24 (6 x 4) lattices. Figs. 1,2. 
The lattice A^ = 18, belonging to the class of lattices A^ = 10, 20, 26, is expected to 
suppress more symmetrical fluctuations, including the four-sublattice state fluc- 
tuations (because A^/2 = 9 is odd), whereas the 6x4 lattice is rather closer to the 



4x4 lattice. It is seen that the tendency is conserved: For the A^ = 18 site lattice 
the self-consistent theory practically reproduces the exact data up to a ^ 0.4, 
whereas for the N = 24 lattice it evidently underestimates the role of frustrations 
starting from a = 0. On the other hand, Hirsch- Tang's consideration [14], which 
does not take into account quantum fluctuations, predicts less sublattice magne- 
tizations (Fig.l) in both kind of lattices. Notice that the latter theory leads to 
completely wrong correlators (they are not presented in Fig.2) for the discussed 
lattices. Therefore, we suggest that the symmetrical lattices A^ = 16, 24 reproduce 
in a more adequate way the main properties of the thermodynamic limit. 

In what follows we address the most symmetrical 4x4 lattice. The scal- 
ing parameter U = f / g (see Ref.4) appears in the self-consistent spin-wave the- 
ory through a Hartree-Fock decoupling of the quartic terms in the Hamiltonian. 
Within the theory, U is given by the self-consistent equations. In principle, one 
can use ?7 as a variational parameter in the spin-wave ansatz 

i,s ~ Pexp(^Va+6tk) \Neel) . (2) 

k 

Here \Neel) is the classical Neel state. The weight factors Wk are defined by t^k = 
V]i/u]f^ , fk and -Uk being the well-known Bogoliubov coefficients; P is a projection 
operator, and the prime means that the sum runs over the small Brillouin zone. 
This variational state is studied in Ref. [19]. Here we treat f/ as a fitting parameter 
obtained from a requirement for best fitting between the sublattice magnetization 
Mg, as obtained from the theory, and the exact-diagonalization results. The 
reasoning for such a consideration comes from the following observations: 

(a) From Takahashi's condition S^ = 0, when applied in the limit A^ = oo, 
one can directly deduce the following scaling relation connecting the sublattice 
magnetizations in the linear spin- wave approximation mo(Q;) (when U = 1) and 
in the self-consistent theory m{a): 

'm{a) = mo{Q:U), N = oo. (3) 

It is interesting to notice that the above scaling relation does not explicitly depend 
on the site spin S (apart from a trivial linear term ). The implicit dependence is 



hidden in the scahng factor U which for a = reads: 

^^1-0.102/25 + 0|(25n ^^^^ 
l + 0.l58/2S + O[{2S)-^] 
For 5* = 1/2 one gets U = 0.775. The next-order term shghtly diminishes the 

latter number. The self-consistent theory predicts a monotonic desrease of U 

versus a in the whole range where the classical Neel state is stable {U ~ 0.6 at 

the phase transition point ac ~ 0.62). If one takes the self-consistent theory as 

a starting point for a systematic perturbation expansion, one can hardly expect 

any drastic qualitative change in the behavior of U versus a. The arguments 

are as follows: First, in this approximation the unphysical modes due to the 

degeneracy of the classical ground state at a = 1/2 acquire gaps, so that some 

of the problems concerning the standard spin wave approach here are resolved. 

Second, the denominator in Eq.(4) is just the rescaling factor of the spin-wave 

velocity, which is expected (also from other methods) to be slightly a-dependent 

and finite at the phase boundary. 

(b) A direct calculation of the scaling factor U ior N = oo and A^ = 16 gives 
practically the same function U{a). In other words, the A^ = 16 lattices is large 
enough in respect to this quantity, as it should be expected, because the factor U 
is a ratio of two short-ranged bosonic correlators. This observation will be used to 
get information concerning the N = oo system. The results coming from an exact 
fitting of the theoretical and exact-diagonalization functions M^ are as follows: 

(i) The resulting scaling factor U is approximatetly a-independent with a 
value close to the one predicted by Eq.(4). (ii) It is seen a remarkable fit to 
the exact correlators (Figs. 3) practically in the whole interval up to « ^ 0.45. 
The misfit for a > 0.45 is easily indicated in the ground-state energy because the 
small overestimates for the short-range correlators, noticed in Fig. 3a, are summed. 
Nevertheless, we have checked that the energy is approximately unchanged by the 
fitting up to q; ~ 0.45. 

Based on the argument (b) and the suggestion that the A^ = 16 lattice better 
refiects the N = oo limit (as compared to less symmetrical lattices), one can 
predict the same picture in the thermodynamic limit: Namely, the monotonic 
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decrease in the scaling factor U vs a should become smoother (as it is for A^ = 16) 
in a more refined approximation starting from the discussed self-consistent theory. 
In particular, the phase boundary a^ ~ 0.62 should drastically move towards 
smaller a ( because U vs a is approximately constant according to the A^ = 16 
fitting). The last two equations, combined with the hypothesis for a smooth 
decrease of the scaling factor U vs a, give the following lowest limit where the 
Neel phase becomes unstable: a* > 0.49, S = 1/2. As a matter of fact, this 
estimate can be slightly increased if one takes the next-order approximation in 
Eq.(2). 

To summarize, a combined approach relying on a comparison with exact- 
diagonalization data and some intrinsic features of the self-consistent spin-wave 
theory ( the scaling relation (3), the smooth monotonic decrease of U{a), and 
the short-range character of the scaling factor U ) lead to the conclusion that the 
theory, when applied to the extreme-quantum system S = 1/2, overestimates the 
stabilizing role of quantum fluctuations. In addition, the same analysis predicts a 
lowest limit a* = 0.49 where the Neel state is destroyed which is quite larger than 
the previous estimate a ~ 0.4 based on the linear spin-wave theory [1,14] and on 
the A^ = 36 exact-diagonalization results [10]. 

An additional understanding of the features of the self-consistent approxima- 
tion can be obtained from the spin ansatz (2). Here we address the square of 
the sublattice magnetization M^, Fig. 4. The exact M^ is compared to: (i) self- 
consistent theory ; (ii) U = 1, i.e., linear theory; (iii) the ansatz ipi, Eq.(2), 
with U = 1. It is obvious that the function M^ (calculated with the ansatz ipi/2, 
U = 1) strictly follows the form of the exact function Mg{a) in a large region up 
to a = 0.5. The main difference between the theory and ips for U = 1 lies in the 
fact that the variational function does not contain unphysical states. Therefore, 
the increasing misfit in M^ (and in the other correlators in Hirsch- Tang's the- 
ory) is predominantly connected with the enhanced role of the unphysical bosonic 
states in the frustrated system (notice that already the linear spin-wave approxi- 
mation in the pure a = system gives a good estimate for the reduced site spin 



mo = 0.303). 



2.2 Violation of the Marshall-Peierls sign rule 

Recent exact numerical diagonalization studies of small lattices [19] show that the 
ground-state wave function of the S = 1/2 Ji — J2 model violates the Marshall- 
Peierls sign rule for sufficiently large a. Here we present results which are based 
on the spin- wave ansatz (2). Originally, the mentioned rule had been proved for 
bipartite lattices with nearest-neighbor interactions [20]. The latter says that the 
ground-state wave function of 5* = 1/2 Heisenberg antiferromagnet reads 

^ = ^(-ira„|n), a, >0, (5) 

n 

where \n) is an Ising state, p„ being the number of, say, up-spins living on, say, 
A-sublattice. Notice that the proof does not work for a system with antiferromag- 
netic next-nearest-neighbor (diagonal) couplings. As a matter of fact, this rule 
is violated, as mentioned above, in the 4x4 lattice provided the frustration is 
strong enough. 

Firstly, let us rewrite the spin-wave ansatz (2) in the form: 

^1 ~ n [l-wir)S^S^_,,]\Neel) , (6) 

ReA 

where the pairing function u'(r) is defined by 

w{r) = — ^'wkcoskr . (7) 

^^ k 

The vector r in Eqs.(6,7) connects sites from different sublattices. 

From the structure of the ansatz it is clear that the sign rule breaks if, and 

only if, the pairing function w{r) changes its sign for some vector r connecting two 

spins which live on different sublattices. For the 4x4 lattice this is just the vector 

r = X + 2y (and the related by symmetry vectors on the lattice). The pairing 

function w(x + 2y) vs a is presented in Fig. 5. For U = 1, w(x + 2y) changes 

sign at a point practically coinciding with the related N = 00 limit, aM = 0.323 



(this is another indication that this symmetrical lattice covers quite well some 
characteristic features of the infinite system). This characteristic point preceeds 
the instability point a* . These observations were based on the spin-wave ansatz 
(2). The predicted weight of Marshall states (Ising states fulfilling the rule) vs 
a is in agreement with the exact result presented in Fig. 6. Quite surprisingly, a 
recent study of the same problem for ground states with larger total-spin quantum 
numbers, Stotai = 1, 2, and 3, [21] indicates a sharp increase of the weight of non- 
Marshall states near the point a ~ 0.52 which is pretty close to limit 0.49 found 
in the present consideration ( see also Ref. 22). 

3 Concluding remarks 

The analysis presented above was based on a combined approach using exact- 
diagonalization data for small lattices and some intrinsic features of the self- 
consistent spin-wave theory. It was directed towards checking the predictions 
of the latter theory for the extreme-quantum system S* = 1/2. It was found a 
stable tendency, namely, the theory excellently fits to the exact data for less sym- 
metrical clusters (A^ = 10,18,20,26), whereas for A^ = 16 (and also A^ = 24) 
this approach evidently underestimates the effect of the frustrations. At the same 
time, Hirsch- Tang's theory, which does not take into account spin-wave inter- 
actions, systematically overestimates the role of the frustrations for each of the 
mentioned lattices. This tendency probably means that some more symmetrical 
(e.g., four-spin) correlations, which are suppressed in less-symmetrical lattices, 
are not properly taken into account in the self-consistent approach. That is way 
the main conclusions are drawn from the comparison with the A^ = 16 lattice 
which is suggested to reproduce better the properties of the N = oo model. The 
comparison gives an estimate a* > 0.49 for the location of the instability point 
a*, which is higher as compared to some previous results (Refs. 1, 8, 10). The 
analysis also shows that in a more refined approximation, using as a starting point 
the discussed theory, a* should move towards smaller values of a (as a matter of 



fact, the fitting based on the A^ = 16 lattice predicts an approximately constant 
rescaling factor U up to the phase boundary, which would mean, if applied to 
the N = oo system, that the instability of the Neel phase should be close to the 
estimate 0.49. Clearly, one needs some additional arguments in favor of such a 
suggestion ( see, e.g., Ref. 22). Finally, it was shown that the spin-wave projected 
ansatz predicts a finite region ( a < 0.323 in linear spin-wave theory) where the 
Marshall-Peierls sign rule is fulfilled in the frustrated N = oo model. For A^ = 16, 
this result was shown to be in accord with the exact-diagonalization data. 
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Captions of figures 

Fig.l: The square of sublattice magnetization vs a for A^ = 18 and N = 24 

lattices. U = 1 (Use) curve corresponds to Hirsch- Tang's (self-consistent) theory. 

The points are the exact-diagonalization data. 

Fig. 2: Spin-spin correlators < SqSr > (R = nx + my) for the A^ = 18 lattice 

in the self-consistent theory. The respective U = 1 curves, which are not drawn 

here, are very bad. The points are the exact-diagonalization data. 

Fig. 3: The correlators < SqSr > (R = nx + my) for the N = 16 lattice. The 

solid lines are the results of the fitting. The other notations are the same as in 

Fig.l and Fig. 2. 

Fig. 4: The square of sublattice magnetization vs a for the N = 16 lattice. The 

curves 1 and 2 represent the spin-wave results for U = 1 and Use, respectively. 

The curve 3 is calculated with the spin- wave ansatz, Eq.(2), for S = 1/2. The 

points are the exact-diagonalization results. 

Fig. 5: The pairing function i(7(x + 2y), as defined by Eq.(7), vs a, U = 1. 

w(x + 2y) vanishes at aM = 0.323 for U = 1 in the thermodynamic limit A^ = oo. 

Fig. 6: The weight of Marshall states vs a. The curves 1 and 2 correspond to 

U = 1 and Use, and are calculated with ipi, Eq.(7). The curve 3 represents the 

exact-diagonalization results, A^ = 16. 
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